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ABSTRACT 

This is the first of two papers reporting observations and analysis of the 
unusually bright (mb = 14.4), luminous {Mb = —25.5), nearby {z = 0.192) 
narrow-line quasar PHL 1811, focusing on the X-ray properties and the spectral 
energy distribution. Two Chandra observations reveal a weak X-ray source with 
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a steep spectrum. Variability by a factor of 4 between the two observations 
separated by 12 days suggest that the X-rays are not scattered emission. The 
XMM-Newton spectra are modelled in the 0.3-5 keV band by a steep power law 
with r = 2.3 ±0.1, and the upper limit on intrinsic absorption is 8.7 x lO^'^cm"^. 
The spectral slopes are consistent with power law indices commonly observed 
in NLSls, and it appears that we observe the central engine X-rays directly. 
Including two recent Swift ToO snapshots, a factor of ~ 5 variability was observed 
among the five X-ray observations reported here. In contrast, the UV photometry 
obtained by the XMM-Newton OM and Swift UVOT, and the HST spectrum 
reveal no significant UV variability. The inferred from the Chandra and 
contemporaneous i/STspectrum is — 2.3±0.1, significantly steeper than observed 
from other quasars of the same optical luminosity. The steep, canonical X-ray 
spectra, lack of absorption, and significant X-ray variability lead us to conclude 
that PHL 1811 is intrinsically X-ray weak. We also discuss an accretion disk 
model, and the host galaxy of PHL 1811. 

Subject headings: quasars: emission lines — quasars: individual (PHL 1811) — 
X-rays: galaxies 



1. Introduction 

The standard model for Active Galactic Nuclei (AGN) and Quasi-Stellar Objects (QSOs) 
proposes that the broad-band optical and UV continuum originates in an accretion disk. The 
X-ray emission is a separate component, produced in a corona located in the vicinity of the 
disk, that creates the X-rays by inverse-Compton-scattering the disk photons. This broad- 
band continuum is then thought to illuminate the gas that forms the broad-line region, 
causing it to emit lines via photoionization. 

Despite the commonality of the origin of the optical through X-ray continuum emission 
and emission lines, there are theoretical reasons that the spectra should vary among indi- 
vidual objects. The origin of differences may be extrinsic; for example, the brightness and, 
to some extent, the shape of the continuu m spectrum of the accretion disk should vary with 



viewing angle (e.g.. iLaor &: Netzer Ill989l ). The origin may also be intrinsic. Even for very 



^ Based on observations obtained with XMM-Newton, an ESA science mission with instruments and 
contributions directly funded by ESA Member States and NASA 
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simple disk models, in which the spectrum is constructed from a sum of annuli locally emit- 
ting as black bodies, the continuum spect rum from the accretion disk should depend on the 
black hole mass and accretion rate 



le.g. 



Frank. King &: Raine Ill992l ). More sophisticated 



accretion disk models that include a range of physical processes ex pected to be important 
also predict a range of shapes (e.g.. lKato. Fukue. fc Mineshige Ill998l ). Another complication 
is that the type of accretion d isk present is predi cted to depend on the accretion rate relative 
to the Eddington value fe.g.. IChen et al. 1119951). and the type of accretion disk can be dif- 
ferent at different radii (e.g.. ISvensson fc Zdziarski Ill994l ). The X-ray emitting corona adds 
another dimension of complication, as its origin and geometry are not very well understood. 
Thus, in principle, the coronal emission may be important or not depending on how much 
of the accretion energy is fu nneled to it, and indeed , we see evidence for a range of coronal 
activity in X-ray novae (e.g.. 



Kubota fc Donell2004J l 



In addition to the theoretical expectation of a range of predicted spectral energy dis- 
tributions among AGN, there is observational evidence that such a range exists. In a study 
of the mult i wavel ength properties of an X-ray-selected heterogeneous sample of quasars, 
Elvis et al. I (1l994j ) observ ed a wide range of spect r al energy distributions fSED s). IWilkes et al 



1994 ). and more recently ISechtold et al~l J2003h , Istrateva et aTl J2005h , and lSteffen et al. 
(l2006l ). found that aox, the point-to-point slope between the 2500A and 2 keV, is inversely 



correlated with the UV luminosity. While some of the large range of spectral energy dis- 
tributions can be accounted for by extrinsic effects such as reddening and absorption (e.g., 
Brandt. Laor. fc Wills 1120001 ). that certainly cannot account for the behavior of all objects. 



PHL 1811 is a nearby {z = 0.192), luminous {Mb = —25.5) narrow-line quasar. 
PHL 1811 was first ca taloged as a blue object in the Palomar-Haro-Luyten plate survey 

( iHaro fc Luyten Ill962l ) . It then rediscovered in the optic al foUowup of the VLA Faint Images 

of the Radio Sky at Twenty Centimeters (FIRST) survey (IWhite et al. 1119971 : iBecker. White, fc Helfand 
19951 ). It is extremely bright (B=14.4, R=14.1); it is the second brightest quasar at z > 0.1 
after 3C 273. Being so bright, it is a very good background source for studies of the in- 
tergalactic and interstellar medium; furthermore, a FUSE obse rvation found i t s spec trum 



to have a rare Lyrn an-limit system that has been studied by iJenkins et al. I (120031 ) and 



Jenkins et al. I (120051 ). It was odd, however, that such a bright quasar was not detected in 
the ROSAT All Sky Survey (RASS). In comparison with other quasars of its luminosity, the 
expected RASS couri t rate is about 0.5 s^^; we placed an upper limit of 1.3 x 10^^ counts s~^ 
(ILeighly et al. Il200ll ). A pointed BeppoSAX observation detected the object, but it was 
still anomalously weak. Too few photons were obtained i n the BeppoSAX ob servation to 
unambiguously determine the cause of the X-ray weakness; iLeighly et al. I (120011 ) speculated 
that either it is intrinsically X-ray weak, or it is a nearby broad-absorption line quasar and 
the X-ray emission is absorbed, or it is highly variable, and we caught it both times in a low 



-4- 



state. 



In this paper and the companion paper ( Leighly et aL 2006|, hereafter Paper II) we 



report the results of several UV and X-ray observations of PHL 1811 designed to explore 
the origin and consequences of the X-ray weakness of this object. First, coordinated Chan- 
dra and ilfST observations were made in 2001. In 2004, an XMM-Newton observation was 
made, and most recently, PHL 1811 was the target of two Swift target-of-opportunity ob- 
servations. In §2 we describe the results and analysis of the Chandra, XMM-Newton, and 
Swift observations, as well as the results of a three-day optical photometry run at MDM 
Observatory. We also compare the XMM-Newton spectrum of PHL 1811 with those from 
other NLSls. In §3 we comment on the long time-scale X-ray and UV variability. We present 
an updated spectral energy distribution in §4. In §5, we discuss the nature of the intrinsic 
X-ray weakness and present an accretion disk model for PHL 1811. We also comment o n 
the apparent spiral host galaxy discovered in the image presented by iJenkins et al. I (120051 ). 
We summarize our findings in §6. Paper II describes the HST and ground-based optical and 
UV observations, and presents Cloudy models that exp lore the unusual emission-line prop - 
erties. Some of the results we r e pres ented previouslv in iLeighly. Halpern &: Jenkins I (12004J ). 
Choi. Leighly fc Matsumoto I (120051 ) , and iPrescott I (120061 ) . We assume a fiat Universe with 
Hq = 70kms~^Mpc~^ and fl^ac = 0.73 unless otherwise specified. 



2. Observations and Analysis 

2.1. Chandra Observations and Analysis 

The Chandra observations were made in imaging spectroscopy mode with the image of 
PHL 1811 placed on the ACIS-S3 detector. The observing log is given in Table 1. We verify 
that the position of the X-ray source is consistent with that of the quasar (Fig. 1). 

The level 2 events files were recreated using the standard procedure. The small correc- 
tion for the time-dependent gain was applied using the corr_tgain program, and the correc- 
tion for the time-dependent ancillary matrix was made using the IDL program acisabs.prc|§. 
The total count rates observed within a circular region 3.94 arc seconds in radius were 
9.7 X 10"'^ counts s~^ and 4.0 x 10~^ counts s~^ from the first and second observations, re- 
spectively. These rates are low enough that pile-up is negligible. Between 0.3 and 9.8 keV, a 
total of 81 and 374 photons were obtained. Based on the background collected from source- 
free areas of the chips, we expect 1 and 3 of these photons to originate in the background. 



^http://www. astro. psu.edu/users/chartas/xcontdir/xcont. html 
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Thus we can conclude that we observe a significant change in fiux in the object by about a 
factor of four between the two observations separated by 12 days. 

A sufficient number of photons were collected in the second observation to look for short 
time scale variability. We used the Bayesian Blocks program available in the ISIS softwar^, 
but did not find any indication of variability during the observation at significance levels 
greater than about 1 sigma. We also tried binning the light curve with 100-second bins, and 
then grouping together by hand points that appeared low or high. The for a constant 
model was 13.5 for 6 degrees of freedom, indicating that variability was marginally detected 
with confidence of 96.4%, although the probability that this variability is significant may be 
lower because the results are biased by preselecting the bin sizes. Note that the A^^ = 6.63 
uncertainty on the fitted constant model is 14%, indicating that we are only sensitive to 
variations larger than this value at the 99% confidence le vel. Such high a mplitude variations 



are rare but not unprecedented in luminous NLSls (e.g.. iLeighly Ill999al ). 



The spectra were accumulated and grouped so there were ~ 20 photons in each bin. 
We first fitted each spectrum between 0.3 and 5 keV separately with a power law and fixed 
Galactic column of 3.76 x lO^'' cm~^, obtained using the HEASARC riu too fl We note that 
the G alactic Lya line in the medium-resolution UV spectrum of PHL 1811 (jjenkins et al. 



20051 ) is consistent with this value of Nh- This model fits both spectra well, yielding photon 



indices of 2.01+°J^ and 2.58lQ;}g for the first and second Chandra observations, res pectively. 



These photon indices are consistent with those observed from NLSls by ASCA (ILeighly 



1999bl ). Note that the ASCA photon indices were taken from models spanning the ~ 0.5-10 
keV band that include a soft excess, warm absorber, and iron line as necessary. Thus, the 
Chandra spectra from PHL 1811, fit over 0.3-5.0 keV (0.36-5.96 rest frame), are consistent 
with the hard X-ray power law found in ASCA spectra of NLSls. 

The best-fitting photon index is steeper for the second, brighter observation, suggest- 
ing that shape of the spectrum has changed between the two observations. However, the 
uncertainties indicate that the difference is not statistically significant. Fitting the spectra 
simultaneously and using the F-test shows that the improvement in the fit represented by 
a change in the photon index is significant at only the 68% confidence level. Fig. 2 shows 
these model fits. 

In order to see if there is any evidence for intrinsic absorption, we next fit the spectrum 
from the second observation with a model consisting of a power law, absorption fixed at the 



ht tp : / /space . mit . edu / CXC / analysis / SITAR/index . html 
'^littp:/ /heasarc. gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nli. pi 
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Galactic value, and absorption in the rest frame of the quasar. We find no improvement in 
the fit and the best-fitting value of additional intrinsic absorption is zero. The Ax^ = 2.71 
upper limit on the intrinsic absorption column is 1.8 x 10^°cm~^. The upper limit is rather 
low despite the poor photon statistics in the spectrum because the spectrum is convex (Fig. 
2). 

The convex residuals of the power law fit to the second observation suggest the presence 
of a soft excess (Fig. 2). Soft excess components are common in the spectra of NL Sls, and 



in the case of poor statististics can be fit adequately by a black body model (e.g., iLeighly 



1999bl ). We add a black body component to the power law model for the second spectrum and 
find that the fit improves by Ax^ = 4.9, and the residuals are fiat (Table 2). However, the 
improvement in the fit is not statistically significant; the F test shows that the improvement 
in fit is significant at only the 71% confidence level. Thus, we cannot conclude that a soft 
excess is present because of the poor statistics. 

The power law index for the power law plus black body fit to the second observation 
is flatter than for the power law alone (2.22 ± 0.34 versus 2.58lo'J8), ^'^^ completely 

consistent with that of the first observation (2.01]'^Q;3g). This suggests that the spectral 
variability originates as an emergence of the black body component when the object is 
brighter. To investigate this possibility, we fit both spectra simultaneously with a power law 
plus blackbody model, fixing the normalization of the black body for the first observation to 
zero. The model fits the data adequately {xl = 0-67 for 16 degrees of freedom (d.o.f.)). The 
jointly-fit photon index is 2.12 ± 0.25, and for this parameterization, the normalizations of 
the power law differ by a factor of 3.5. 

Spectral ratios provide a complementary and model-independent approach to the ques- 
tion of spectral variability. We rebin the spectrum from the second observation to the binning 
of the first spectrum and take their ratio (Fig. 2 inset). This shows that the spectral vari- 
ability is predominately in the softest band, supporting our hypothesis that a variable soft 
excess is responsible for the spectral variability. In this case, the variability of the power 
law component is a factor of 3.5. The ^ constant ratio model is 3.8 for 3 degrees of 

freedom which significant at the 71% confidence level. 

To summarize the results of the Chandra observations, we find conclusive evidence for 
factor of ~ 4 variability between the two observations separated by 12 days. Detailed anal- 
ysis is hampered by poor statistics; however, we find marginal evidence {2a) for variability 
on times scales of thousands of seconds in the second observation when the object was 
brighter. Spectral fitting reveals a steep spectrum with no evidence for intrinsic absorption, 
and marginal evidence for spectral variability between the observations, with the spectrum 
becoming steeper when the object is brighter. The spectral fitting and the ratio of the 
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spectra suggest that the spectral variabihty is caused by the emergence of a soft excess com- 
ponent when the object is bright. The measured spectral indices range between 2.0 and 2.6, 
depending on the mod el. They are con sistent with the photon indices measured in ASCA 
observations of NLSls (ILeighly Ill999bl ). This fact suggests that we see the intrinsic X-ray 
emission from the central engine, and that the X-rays are powered by inverse Compton 
scattering of soft photons as in other AGN. 



2.2. XMM-Newton Observation and Analysis 



PHL 1811 wa s observed by XMM- Newton 1 November 2004 using the EPIC PN flStriider et al. 



2001 ) and MOS (ITurner et al. Il200ll ) instruments and the optical monitor (IMason et al. 
200ll ). The EPIC observations were carried out using the Thin filter in PrimeFullWindow 
mode. The data were reduced using standard selection criteria. The object was observed for 
32.1 ks using the MOS detectors and for 27.5 ks using the PN. The details of the observation 
are given in Table 1. 

Background flares are a concern in XMM-Newton data analysis, and we observed the 
background to vary during the observation. For most of the observation, the background 
relatively low and stable. However, even at the lowest rate, it appears to be slightly elevated 
compared with the quiescent rat^ by a factor of approximately 2.3 for the PN and 1.4 in the 
M0S1+M0S2 in both the soft (0.5-2.0 keV) and hard (2.0-10 keV) bands. In addition, for 
the first 3000 seconds in the PN detector, and the first 5000 seconds for the MOS detector, 
there occurred a small background flare that was higher than the quiescent rate by a factor of 
2-3. We extracted and analyzed spectra with and without this flaring period, and conclude 
that the flare is so small that it does not adversely affect the results. Therefore, we analyze 
the entire exposure. 

We extract light curves from the PN and MOS detectors. The target is relatively 
weak, so we use a source extraction region with a radius corresponding to an encircled 
energy function of 80%; the radius was 27" for the PN and 23" for the MOSl and M0S2. 
Background light curves and spectra were extracted from nearby, source-free regions of the 
detector. The background spectrum is flatter than the source spectrum, and we find that the 
background contribution to the emission in the extraction aperture is equal to that of the 
target at ~ 5 keV. The background dominates at higher energies; therefore, we extract light 
curves in the 0.5-5 keV band. We use a bin size of 1000 seconds, which yields an average of 
~ 30 source counts per bin, and the background contributes about 12% of the total counts. 



^XMM-Newton User's Handbook §3.3.8 |http:/ /heasarc. gsfc.nasa.gov/docs/xmm/uhb/node38. html 
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The resulting background-subtracted light curve is shown in Fig. 3. The mean count rate 
is 2.8 ± 0.17 X 10"'^ counts s"^. The light curve is statistically consistent with no variability 
{x^ = 33.4 for 30 d.o.f. for a constant model). Note that the Ax^ = 6.63 uncertainty on 
the fitted constant model is 9%, indicating that we are only sensitive to variations larger 
than this value at the 99% confidence le vel. Such high a mplitude variations are rare but not 
unprecedented in luminous NLSls (e.g.. iLeighly Ill999al ). 



The spectra were extracted using the regions described above, and grouped so that there 
were 15 photons per energy bin. We fit the PN in the 0.3-5 keV range and MOS in the 0.5-5 
keV band, simultaneously. The spectra yield 611, 191 and 204 source photons for the PN, 
MOSl and M0S2. The spectra are fit very well with a model consisting of a constant, a 
power law and Galactic absorption column described in §2.1. The results are listed in Table 
3, and the model fit is shown in Fig. 4. Note that unlike the second Chandra observation, 
the residuals are flat, and there is no evidence for a soft excess component. 

The spectra are adequately fit using a steep po wer law (F = 2.3 ± 0.1), typical of 
that from NLSls observed by ASCA (ILeighly Ill999bl ). There is no evidence for additional 
absorption. We add a neutral absorption component at the redshift of PHL 1811 to the 
model but find no significant reduction in x^- The 90% confidence upper limit (Ax^ = 2.7) 
on additional absorption is Nh = 8.7 x 10^°cm~^. 

The constant in the model is fixed to a value of 1 for the PN spectrum, and allowed 
to be free for the MOSl and M0S2 spectra. It can be seen in Table 3 that the best- 
fitting normalizations of MOSl and M0S2 spectra are 25% and 32% higher than that of 
the PN. This difference cannot be explained by residual cali bration uncertainty bet ween 
the XMM-Newton EPIC instruments as that is now quite low (jStuhlinger et al. II2006I ). We 
analyze the X-ray spectra from another object in the field (located at RA=21h54m41s, 
Dec=— 9d26h49m) that is about 50% brighter than PHL 1811. The normalizations of the 
EPIC spectra for this object were completely consistent with one another (MOSl constant: 
1.00±0.12; M0S2 constant: l.OOl^J;}^). Another object in the field of view was independently 
analyzed by another of the authors (D. Grupe) with the same result. 

While the small numbers of photons in the PHL 1811 spectra mean that the difference 
in normalizations among the models for the different EPIC spectra is not statistically signif- 
icant, there still seems to be a problem with the spectra that needs to be understood. We 
believe that the problem originates in optical loading. The nom inal lirnit for optical loadiri g 
using the thin filter is V ?5i 12 for both the PN and the MOS JSmith 1 12004 : lAltieri |[2003[ l. 
PHL 1811 is a fainter optical source (B=14.4, R=14.1); however, it has a very blue spectrum 
and is a very weak X-ray source, and thus is an unusual object compared with the stellar 
calibration sources used to determine the loading limits. Another point that supports our 
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contention that optical loading is important is the fact the optical loading in the M0S2 is 
observed to b e about 7% la rger than that in the MOSl, possibly due to variations in the filter 



transmission (lAltieri II2003I ): we also observe a higher normalization for the M0S2 spectrum. 
However, the degree of contamination by optical photons cannot be very large because the 
MOS spectra do not show any observable distortion; a power law fit to them alone yields 
an identical photon index as that obtained from the PN. Nevertheless, given the fact that 
the MOS spectra are possibly contaminated by optical loading, we measure the flux from 
the PN spectrum, noting that we cannot be sure that this spectrum is uncontaminated by 
optical loading as well, and may represent an upper limit on the flux for this observation. 

PHL 1811 was observed using the Optical Monitor with the UVM2 filter. Ten exposures, 
each with duration of 2580 seconds, were made. The SAS task omichain was run to reprocess 
the data. The count rate information was extracted from the omichain output. In addition, 
the count rates were extracted from the images using the IRAF task phot following the proce- 



dure described on http: / /xmm. vilspa.esa.es/sas/documentation/watchout/uvflux.shtml The 



count rates were then converted to flux using the conversion factor for the UVM2 filter 
(2.17 X 10~^^ erg s~^ cm~^ A~^ count"^). The results from the two extraction procedures were 
consistent to within 2%, and therefore hereafter we discuss the results from the omichain 
output. 

One of the goals of the OM observation was to look for UV variability. A constant model 
fit to the light curve yielded = 16.86 for 9 degrees of freedom. Thus, a constant model is re- 
jected at the 94.9% confidence level. However, we do not consider this evidence for marginally 
significant UV variability because the fluctuations are different for the omichain and IRAF- 
reduced data. The mean flux of the 10 observations was 2.773±0.007x 10"^"^ erg s~^ cm^^ A^^. 



2.3. Joint Chandra and XMM-Newton Modelling 

In principle, better constraints on spectral variability can be obtained by jointly fit- 
ting the Chandra and XMM-Newton spectra. We first fit with a power law plus Galactic 
absorption model, allowing the normalizations to be free among the two Chandra spectra 
and the XMM-Newton PN spectrum, but with the photon indices constrained to be equal. 
The spectra are fit adequately; the is 70.5 for 70 degrees of freed om. The best fitting 



photon index is 2.36;!;q :[^, again typical of an NLSl (iLeighly I Il999bl ). We notice that the 
normalization for the spectrum from the first Chandra observation is identical to that of the 
XMM-Newton PN spectrum; if we tie these two parameters together, there is no change in 
(x^ = 70.5 for 71 degrees of freedom). These two spectra seem to describe the object 
in the same state (see also §3), so we leave their parameters tied together throughout the 
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remainder of this section. 

Next, we allow the photon index of the second, higher flux, Chandra observation to 
be flt independent of that of the flrst Chandra spectrum and the XMM-Newton spectrum. 
We obtain a somewhat better flt with = 64.1 for 70 degrees of freedom. However, the 
decrease in is signiflcant at only the 63% level according to the F test. The resulting 
photon indices are 2.22 ±0.14 for the low state (flrst C/ian(ira observation and XMM-Newton 
PN spectrum), and 2.57^0^8 for the second, brighter Chandra observ ation. These in dices 



are still within the range of power law indices observed from NLSls (jLeighly Ill999bl ). In 



addition , it is of ten found the NLSl spectra soften as when they are brighter (e.g.. Fig. 4 of 



Leighly I (jl999a] )): thus, intrinsic softening of the spectrum is a plausible explanation for the 



marginal spectral variability we observe. 

It is possible that the variability results from variable cold absorption. We test this 
scenario by constraining the photon indices and normalizations to be the same for all three 
spectra, and including a neutral absorption column in the quasar rest frame in the model, 
allowing the absorption column to vary. This model gives a very poor flt; the reduced 
= 2.57 for 70 degrees of freedom. Thus, we can reject on statistical grounds the idea that 
the spectral and flux variability originates solely in variable neutral absorption. Allowing the 
normalizations to also be free yields a good flt (x^ = 65.3 for 69 degrees of freedom), although 
it is not a signiflcant improvement over the no-absorption model ( Ax^ = 5.2, signiflcant at the 
58% confldence level according to the F test). The additional absorption is consistent with 
zero for the brighter Chandra observation with an upper limit of 1.4 x 10^° cm~^, and is equal 
to5.0l|^xl020 cm for the fainter Chandra observation and the XMM-Newton observation. 
Although this model (variable cold absorption plus variable power-law normalization) flts 
the spectra well, the scenario seems unlikely because it requires that absorber variability be 
coordinated with intrinsic flux variability. 



2.4. The X-ray Spectra of NLSls; How does PHL 1811 Compare? 



ASCA observations of NLSls demonstrate that their spectra can generally be described 
by a hard power law with average index of 2.19 ±0.10 and a soft excess that can be modeled 
using a blackbody component. The strength of the soft excess varies from object to object, 
with the o bjects haying t he overall steepest spectra also showing the highest amplitude 
variabihty (jheighly 1999b ). 



XMM-Newton, with its large effective area, has revolutionized our understanding of the 
X-ray spectra of AGNs. Soft excesses are now seen to be relatively common in quasars 
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in general (e.g., iPorquet et al. 1 12004| ). NLSl spectra are still found to have soft excess 
that sometimes can be modelled by a dual Comptonization model, in which soft photons 
are scattered by Comptonizing media of two different tempera tures. Example s of objects 



with this type of spectrum are Ton S180 (jVaughan et al. II2002I ) and Mrk 896 (jPage et al. 



20031 ). In other cases, the X-ray spectrum is very complex, with a very prominent soft ex- 
cess and complex abs orption features at h igh energies. T hese spectra can be mode lled using 
partial covering (e.g. iTanaka et al. 1 12004| ) or refl ection (IMiniutti fc Fabian I l2004r). Exarn - 



ples of this type of object include IH 0707— 495 (IGallo et al. 



IRAS 13224-3809 flBoUer et al. II2003I ). and PHL 1092 fIGallo et al. 



2004 



3 



Tanaka et al. I l2004f) 



20M) 



Gallo 



(120061 ) 



split NLSls into two classes: those with and without significant complex features in their 
high energy spectra. He finds that objects with complex high energy spectra tend to be X- 
ray weak, and proposes that the X-ray weakness is consistent with either attenuation in the 
partial covering scenario, or the focusing X-rays away from our line of sight in the reflection 
scenario. 

How does the XMM-Newton spectrum of PHL 1811 compare with those from other 
NLSls? It is important to note that the quality of our spectrum is much poorer than those 
from many XMM-Newton observations of NLSls, such as those mentioned above, due to its 
low X-ray flux. The PN spectrum has 857 photons between 0.3 and 5 keV in the observed 
frame, corresponding to 0.36-5.96 keV in the rest frame. Of these, 246 are background, 
leaving 611 net source photons. This spectrum is adequately modelled with a power law 
and absorption originating in our Galaxy. The reduced is 0.94 for 51 degrees of freedom, 
and the photon index is 2.25 ± 0.15. The ratio of the data to the power law plus Galactic 
absorption model is shown in Fig. 5. This figure shows that there are no residuals that 
suggest a more complex model, and indeed, because the reduced is less than one, a more 
complex model would over-parameterize the data and would not be justified statistically. It 
is important to note that this photon index is consistent w ith the average hard X-ray photon 
index observed in ASCA spectra from NLSls (2.19 ifa 0.10: iLeighlv lE999b[ ). Thus, PHL 1811 
resembles an average NLSl without a soft excess. The power law X-ray spectrum is believed 
to be produced by Compton upscattering of soft photons in a hot plasma. This is the same 
process believed to operate in AGN in general; the reason that the photon index is steeper 
in NLSls than in broad-liiie qua sars because the hot plasma has been Compton cooled (e.g.. 
Pounds. Done fc Osborne Ill995l ). 



Is the lack of complexity in the PHL 1811 spectrum due to the low flux and poor 
statistics, or does it really have a simple spectrum? We can obtain some answers to this 
question by extracting sufficiently short segments of XMM-Newton data from other NLSls so 
that we obtain spectra with approximately 611 photons between 0.3 and 5 keV. We perform 
this exercise on two NLSls: Ton S180 and IH 0707-495. 
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Ton S180 has a spectrum with a mild soft excess. IVaughan et al. I (120021 ) model it 
using a dual Comptonization model, and it was classified by iGallo I (120061 ) as an NLSl with 
a "simple" spectrum. The Ton S180 data are seen in the middle panel of Fig. 5. This 
object is so bright that we need an exposure of just 51.52 seconds to obtain 572 photons 
between 0.34-5.61 keV (the observed-frame range corresponding to the rest-frame range of 
0.36-5.96 keV for this z = 0.062 object). We fit the spectrum with a power law plus Galactic 
absorption, and obtain a good fit, with a photon index of 3.041° ;}g and xl = 0-99 for 33 
degrees of freedom. The fact that the reduced is less than one means that there is no 
statistical evidence for spectral complexity; in addition, we see no suggestive residuals in the 
data-to-model ratio. However, in contrast to PHL 1811, the photon index is significantly 
steeper than the average from NLSls observed by ASCA, clearly because in this spectrum 
we are fitting the soft excess component predominantly. If this were the only X-ray spectrum 
that we had of this object, we would suspect that spectr al complexity may be present and a 
hard tail might be seen in a longer observation, as it is (IVaughan et al. 1120021 ). 



IH 0707—495 has a complex spe ctrum that wa s modeled by iGallo et al. I (l2004bl ) using 
partial covering, and was classifie d bvlGallo I (20061 ) as an NLSl with a "complex" spectrum. 
Note that this is the class that iGallo I J2006f r observe to be somewhat X-ray weak. The 
IH 0707—495 data are seen in the lower panel of Fig. 5. In this gment 289 seconds 

in length yielded a spectrum with 602 photons between 0.34 and 5.73 keV (the observed- 
frame range corresponding to the rest-frame range of 0.36-5.96 keV for this z = 0.041 object). 
We fit the spectrum with a power law plus Galactic absorption. In this case, the resulting 
photon index is even steeper than for Ton S180 (F = 4.17!^q]^7), reduced is significantly 
greater than 1 (x^ = 1.61 for 34 degrees of freedom), and significant residuals are seen in the 
ratio. Clearly, despite the poor statistics in the spectrum, the need for a complex model is 
evident. It is also clear that although this object is X-ray weak, the X-ray spectrum is very 
dissimilar to that of PHL 1811. 

This exercise shows that the X-ray spectrum of PHL 1811 is clearly different than 
those of two other NLSls, Ton S180 and IH 0707—495. These objects are representative 
of the complexity of spectra from NLSls observed by observed by XMM-Newton. The 
XMM-Newton spectrum from PHL 1811, a simple power law with a photon index consistent 
with the mean hard X-ray index from NLSls observed by ASCA, suggests that the X-ray 
emission mechanism is simple Compton-upscattering of soft photons in a hot plasma, as 
in other AGN, and that the spectrum is unaltered by any extrinsic effects such as partial 
covering or reflection. 

It may also be important that PHL 1811 is more UV-luminous than other NLSls. The 
monochromatic luminosity at 2500A is 30.9; (Paper II); thus it is seen to be about 5 times 
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more luminous than the most lumi nous obje c t show n i n Fig. 2 of iGallo I (120061). Also, note 



that a very similar figure is shown in iLeighly I (120011 ) and lMatsumoto. Leighly fc Kawaguchi 



2.5. Swift Observations and Analysis 



PHL 1811 was observed by the Swift Gamma- Ray Burst Explorer Mission (iGehrels et al 



20041 ) on 2005 October 22 starting at 09:52 UT for a total of 2.5 ks. The observ ations were 
performed simultaneously with the X-ray Telescope (XRT; iBurrows et al. Il2005r) in t he 0.3- 
10.0 keV energy range and the UV-Optical Telescope (UVOT; iRoming et al. II2005I ) in the 
1700-6500A wavelength range. It was observed again on 2006 May 12 for 1.6 ks. The details 
of the observations are given in Table 1. 

The XRT data reduction was performed by the task xrtpipeline version 0.9.9, which 
is included in the HEAsoft package 6.0.4. Source photons were selected in a circle with a 
radius of 23.4" and the background photons in a source-free region close by with a radius of 
95". Those photons were extracted and read into separate event files with XSELECT version 
2.3. Twenty-two photons were detected during the first observation, and ten were detected 
during the second observation for count rates of8.8±1.9xl0 
respectively. 



"^s-i and 6.2 ±2.0 x 10' 



In the 2005 October 22 observation, the UVOT photometry was performed with the 
UV filters UVWl, UVM2, and UVW2. During the 2006 May 12 observation, all six (optical 
and UV) filters were used. The details are given in Table 1. After the aspect correction the 
exposures in each filter were coadded into one image with uvotimsum and the magnitudes 
and fiuxes in each filter were determine with the task uvotsource. The results (observed 
fiuxes, uncorrected for Galactic reddening) are listed in Table 4. 



2.6. MDM Optical Photometry 

Optical photometry data were taken on PHL 1811 at MDM Observatory using the 1.3 
McGraw-Hill telescope on the nights of 2004 October 14, 15 and 1 6 as part of a proi ect to 
search for optical variability in the narrow-line quasar PHL 1092 (IGallo et al. II2004I ). We 
used the thinned CCD "Templeton" and the telescope in the f/7.6 configuration, yielding 
a angular size of 0.50"/pixel. The weather was good on October 14, with typical seeing of 
1.7", although the sky was not photometric. The seeing was worse on October 15 (average of 
2.3") and it was intermittently cloudy. On October 16, the weather had deteriorated further. 
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and few usable frames were obtained. 

We observed PHL 1811 using the I, V, B and U filters, obtaining several exposures 
in each filter. The total number of frames analyzed were 9, 19, 12 and 11 for the I, V, B 
and U filters, respectively. Within each night, the observations were all obtained within a 
time span of about 30 minutes, so we could test the interday optical variability principally 
between October 14 and October 15, as the sampling and data quality was much worse on 
October 16. 

The images were reduced using standard IRAF procedures. Aperture photometry was 
performed on PHL 1811 and 4 field stars using an aperture twice the size of the image PSF 
FWHM. The aperture was chosen to ensure that essentially all of the photons from the 
object were measured, especially in cases where the image was slightly trailed due to clouds 
and loss of tracking. The ratio between PHL 1811 and the field stars was computed and 
errors were propagated. PHL 1811 is a bright object, so the mean signal-to-noise ratios in 
these ratios are high, ranging from 50 to 210. 

We tested for variability in the ratio light curves using the "excess variance" , a technique 
commonly used in X-ray astronomy. The excess variance is a measure of the observed 
variance in the light curve subtracting the variance due to measurement errors. In nearly all 
cases, this measure was negative. We conclude that no evidence for optical variability was 
found. 



3. Long Time-scale X-ray and UV Variability 

PHL 1811 has now been observed in the X-ray bandpass se ven times between 1 990 and 



2006: during the ROSAT AW Sky Survey, by BeppoS AX in 2000 fiLeighlv et al. 1120011 ) . and in 
the two Chandra observation, the XMM-Newton observation, and the two Swift observations 
reported here. In Fig. 6 we show the long-term light curve composed of the secure measure- 
ments, i.e., the last five observations. We do not plot the RASS upper limit or the BeppoS AX 
observation that, with the large detector PSF, was certainly contaminated by X-ray emission 
from neighboring objects. For the Chandra and XMM-Newton observations, the X-ray fiux 
densities were estimated from the best-fitting power-law models, and the uncertainties were 
obtained by propagation of the errors on the normalizations and photon indices. For the 
Swift observations, the fiux densities were estimated from the count rates using PIMM^, 
assuming the XMM-Newton photon index (P = 2.3) and Galactic absorption, and the un- 



ht tp : / /cxc . harvard . edu/ toolkit /pimms . j sp 
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certainties in the flux densities were assumed proportional to the uncertainty in the count 
rate. 

We find that PHL 1811 has varied significantly by a factor of ~ 5 in this time period. 
The first Chandra observation and the XMM-Newton observation found it to be in a relatively 
low state, with vFy at 2 keV rest frame equal to ~ 1.0 x 10~^^ergs~^ cm~^; the fluxes in 
these two observations are consistent with one another. The second Chandra observation 
and the two Swift observations show a significantly higher flux by a factor of 3.5-5.5, and the 
fluxes of these three are roughly consistent with one another, although the uncertainties are 
evaluated differently. Although there are only five points, these data suggest that the X-ray 
flux oscillates between two states that differ by a factor of 4-5. The NLSl IH 0707—495 
seems to behave the same way, oscillating between two flux states that differ by a factor of 
~ 10 jLeighlv et al. Ibood ) 



As will be discussed in J4, quasars with PHL ISll's optical luminosity are statistically 
expected have values of aooLl equal to —1.6. The dashed line in Fig. 6 shows the predicted 
X-ray flux for = —1.6 based on the UV flux of the HST spectrum. Although PHL 1811 
varies significantly, it never approachs the nominal X-ray flux for an object of its UV lumi- 
nosity. 

We have four epochs of UV observations: the HST STIS spectroscopic observation made 
2001 December 3 (discussed in detail in Paper II), the XMM-Newton OM observation made 
2004 November 1, and the two Swift UVOT observations made 2005 October 22 and 2006 
May 12. We search for possible UV variability by comparing the latter three photometry 
measurements with the HST spectrum. 

The effective wavelength of the XMM-Newton OM UVM2 filter is 2310A. We plot the 
observed fluxes on the observed-frame merged HST and optical spectrum from Paper II 
in Fig. 7. We find that the photometry flux is completely consistent with the observed 
spectrum. 

The effective wavelengths of the Swift UVOT photometry points are listed in Table 1, 
and the inferred fluxes are plotted on Fig. 7. Note that the filter transmission functions are 
different in the Swift UVOT compared with the XMM-Newton OM, even though the filter 
names are the same; hence the effective wavelengths are somewhat different. Like the XMM- 
Newton OM photometry, the correspondence between the UVOT photometry and the HST 
spectrum is very good with the exception of the photometry using the UVW2 filter. That 



^aox is defined as the point-to-point slope between 2500A and 2 keV; i.e., aox = log(-F250o/-F2fcev)/2.61. 
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filter is especially difficult to calibrate as the complete ffiter transmission curve is not knowio, 
and because of the paucity of suitable calibration targets. Therefore, we do not consider the 
disagreement with the HST spectrum indicative of a spectral change. We conclude that we 
observe no evidence for any UV variability between the four UV observations that span a 
time period of 4.5 years. 

To summarize, PHL 1811 has now been observed seven times in the X-rays over a period 
spanning 16 years. During this time, the X-ray flux has been observed to vary by a factor of 
~ 5, but it remains well below that of a typical quasar of its UV luminosity. In contrast, the 
four UV observations made over a period of 4.5 years do not show any convincing evidence 
of UV variability. 



4. The Spectral Energy Distribution 



We presented the ffist spectral energy distribution of PHL 1811 in lLeighly et al. I (120011 ) . 
based on an optical spectrum, a i^OS'^T upper limit, a BeppoSAX ohseivaXioTi, and multi- 
wavelength photometry. In Fig. 8, we present an updated spectral energy distribution. For 
the optical and UV, we used the merged spectrum described in Paper II. The Chandra results 
are represented by regions on the SED plot that were constructed using the third joint model 
fit (Table 2). The contours were constructed by successively setting each variable parameter 
to its Ax^ = 2.71 value and computing the model, then determining the maximum and 
minimum of all the models. 

From the best fit Chandra model and the merged HST and optical spectrum, we compute 
to be —2.40 for the first observation, which was nearest in time to the ilfST observation, 
and —2.19 for the second observation. As seen in Figs. 6 and 7, the Chandra and HST fiuxes 
are comparable with the XMM-Newton and Swift fiuxes, so we confine our discussion to the 
Chandra and HST data here without loss of generality. 



Wilkes et al. I (119941 ) compute a regression between optical luminosity and aox for a 



heterogeneous sample of quasars. Using their cosmology {Hq = 50kms~^ Mpc~^, go = 0), 
we obtain a luminosity distance for PHL 1811 of 1261 Mpc, and a corresponding luminosity 
density at 2500A of 1.44 x 10^^ ergs~^ Hz~^. Then, using their regression, we predict aox to 
be —1.6. We plot the predicted X-ray fiux assuming this value of aox on Fig. 8, as well as a 
vertical bar that indicates the range of X-ray luminosities observed by lWilkes et al. I (119941 ). 



More recently, IStrateva et al. I (120051 ) and ISteffen et al. I (120061 ) updated the aox re- 



http: / / swift.gsfc.nasa.gov / docs /heasarc/caldb /swift / docs /uvot / uvot_caldb_filtertransmission_02.pdf 
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gression using a large sample of optically-selected active galaxies that span a large range 
in redshift and luminosity, yet have highly complete X-ray data. For their cosmology 
{Hq = 70kms~^ Mpc"^, Qm = 0.3, and Qa = 0.7) we obtain a luminosity distance of 
936 Mpc, and a corresponding luminosity density at at 2500A of 7.92 x 10^°ergs~^ Hz~^. 
Their regression yields a predicted aox of —1.60 also. From their Fig. 4, we find that the 
envelope of aox observed for quasars of this luminosity spans approximately —1.75 to —1.4. 
Our observed X-ray luminosity densities are factors of 130-450 below the high value, and 
13-45 below the low value. Thus, PHL 1811 is observed to be significantly X-ray weak 
compared with other quasars. 

Brandt. Laor. &: Wills ( 2000 ) c ompil e the distribution of from the PG quasar 



sample studied by lBoroson fc Green I (119921 ). They find a suggestion of a bimodal distribution 
with 10 of the 87 objects classified as X-ray weak with aox < —2. Then, they find there 
to be a connection between aox and the presence of significant C IV absorption lines, such 
that most of the the soft X-ray weak objects have absorption-line equivalent widths greater 
than 5A. They infer these results to imply that X-ray absorption is the primary origin of 
soft X-ray weakness in AGN. 



Clearly, PHL 1811 does not follow the trend observed by lBrandt. Laor. fc Wills I (120001 ). 
It is soft X-ray weak, but, as discussed in Paper H, there is no evide nce for any significant 
intrinsic C IV absorption lines. In the iBrandt. Laor. fc Wills I (120001 ) sample, objects with 
similar aox as observed from PHL 1811 have a C IV absorption line equivalent widths of 
5-20 A. 

Furthermore, the X-ray spectrum shows no evidence for intrinsic absorption. If low- 
column-density absorption were present, we expect to observe a fiat spectrum; in contrast, 
for a single power law model, we measure photon indices of 2-2.6, consistent with unabsorbed 
quasars, and an upper limit on intrinsic absorption of 8.7 x 10^°cm~^. 

Another possibility is that a high-column-density or Compton-thick absorber is present 
in our line of sight, so that we see no direct continuum emission. Our spectra do not probe 
high enough energie s to see whether there i s a highly absorbed c omponent, as has been found 



in BALQSOs (e.g., [Gallagher et al. I l2002l : iGreen et al. Il200ll ). If the continuum emission 



were completely absorbed, then the X-ray spectrum that we see might have been scattered 
into our line of sight (i.e., similar to a Seyfert 2 galaxy or BALQSO). Electron scattering is 
energy independent, so in this scenario we would expect to see the intrinsic power law with 



terandt. Laor. fc Willsl |2000l ) use an alternative definition of Uox using the flux density at 3000A rather 
than 2500A. They estimate that aoa; (2500) — 1.03aoa;(3000) — 0.03a„, where a„ is the slope of the spectrum 
between 2500 and 3000A. 
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attenuated flux. In Seyfert 2s, the electron scattering occurs over an extended area, and 
thus while the intrinsic X-ray emission may vary, the scattered emission does not because 
variability is washed out as it scatters over the extended region. The fact that we see 
significant variability between the two Chandra observations, separated by 12 days, argues 
that we are not seeing scattered light. This conclusion depends on the compactness of the 
electron-scattering mirror; if unusually compact, observation of variability would be possible. 
Regardless, the observation of significant variability between two observations separated by 
twelve days from this luminous quasar suggests that we are seeing the intrinsic emission from 
the AGN. 

Narrow- lin e Seyfert 1 galaxie s are known for their high-amplitude X-ray variability (e.g., 
Leighly Ill999a] ). As discussed in iLeighly et al. I (l200ll ). it was possible, at that time, when 
we had only the RASS upper limit and the BeppoSAX data, that we had coincident ally only 
observed PHL 1811 while it was in a transient low state. In this paper, we report five more 
X-ray observations, all of which find it to be a significantly weak X-ray source. Thus, the 
probability that we coincident ally observe it in a low state is decreasing. PHL 1811 appears 
to be intrinsically X-ray weak. 



5. Discussion 
5.1. PHL 1811 is Intrinsically X-ray Weak 



In iLeighly et al. I (120011 ). we reported the first X-ray detection of PHL 1811 by Bep- 
poSAX. That observation showed that PHL 1811 appeared to be X-ray weak, but the 65- 
net-photon spectrum was not sufficient to determine the origin of the X-ray weakness. We 
presented three alternatives for the X-ray weakness: 1. PHL 1811 is a BALQSO, and the 
X-ray emission is absorbed; 2. since PHL 1811 is an NLSl, it is highly X-ray variable, and 
we happened to catch it in a low state; 3. PHL 1811 is intrinsically X-ray weak. The HST 
observation discussed in Paper II and Chandra observations reported here show that it is 
not a BALQSO, as there is no evidence for UV absorption lines. Furthermore, there is no 
evidence for absorption in the X-ray spectrum, and the significant variability between the 
two Chandra observations suggests that the X-ray emission is not scattered. So the first 
hypothesis is firmly ruled out. 

We can never conclusively rule out the second hypothesis, that PHL 1811 is highly X-ray 
variable and we always just happen to catch it in a low state. However, it has now been 
observed seven times between ~ 1990 (during the ROSAT All Sky Survey) and 2006 May 
(in a Swift observation), and it has never been observed to be bright. In fact, since it has 
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already varied by a factor of ~ 5 among the observations reported here, it may have aheady 
been as bright as it can get. It seems increasingly unlikely that it will ever be as X-ray bright 
as other quasars. 

Thus, we conclude PHL 1811 is intrinsically X-ray faint. Most quasars are bright X-ray 
sources, so wh at property of t he PH L 18 11 central engine ca uses it to be X-ray faint? As we 
pointed out in lLeighly et al. I (120011 ) (also lGrupe et al. Il200ll ). there is no obvious reason why 
intrinsically X-ray weak quasars should not exist. Briefly, a quasar arguably cannot exist 
without accretion as a source of fuel. In an object like PHL 1811, that accretion probably 
occurs through an optically thick, geometrically thin accretion disk that emits the observed 
strong optical and UV continuum (Fig. 8). Such a disk would never be hot enough to emit 
X-rays, which are probably emitted by the corona, a separate component. We postulated 
that there may be situations in which the corona may not exist or may be weak. 



Why might the corona be weak? As di s cussed in iLeighly et al. I (120011 ). PHL 1811 is a 
very luminous NLSl. As discussed by iLaor I (120001 ). narrow Balmer lines in luminous AGNs 
may imply high accretion rates if the width of the lines is dominated by virial motions. 
One possible central engine geometry considers the X-rays to be emitted by a central, hot, 
optically thin, geometrically thick disk, and the optical an d UV emitted by an opticall y 



thick, geometrically thin disk with a large inner radius (e.g.. iZdziarski fc Gierlinski II2004I ). 
At high accretion rates, it might be expected that the inner radius of the optically thick, 
geometrically thin disk would shrink down toward the innermost stable orbit, with the volume 
of the central hot X-ray emitting region shrinking with it, and the intensity of the X-ray 
emission corr espondingly decreasing . This is thought to happen in Galactic X-ray transient 
objects (e.g.. iKubota fc Done 112004 ). 



Alternatively, the corona may lie on top of the optically thick, geometrically thin accre- 
tion disk and be fed by reconnecting magnetic flux tubes buoyantly emerging from the disk. 
Decreasing the amount of energy released by reconnection in PHL 1 811 would result i n wea k 
X-ray emission. Why would that happen? One model, proposed by lBechtold et al. I (120031 ) . 
explains the dependence of on luminosity in the context of a disk/corona model in which 
the two phases are thermally coupled, and in which the amount of energy dissipated into the 
corona depends on the gas pressure in the disk. In this model, is steeper for larger black 
holes, larger accretion rates with respect to Eddington, and larger viscosity parameters. 

It is also possible that since the optical-UV spectrum is very soft, softer than most 
quasars, the corona is flooded by soft photons, catastrophic ally Co r apton cooling it, and 
so reducing the X-ray emission. This idea is discussed by iProga I (120051 ) . who suggests 
that a luminous accretion disk can simultaneously drive an outflow and quench the corona. 
This idea is supported by the fact that we see blueshifted high-ionization lines in PHL 1811 
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(discussed in Paper II), and by the observed inverse c orrelation between aox and the blueshift 
of C IV lines among NLSls Jleighlv fc Moore 112004 ). 

Another way that a high accretion rate may lead to reduced X-ray emission is through 
so-called "photon trapping". If the accreting gas is simply free-falling into the black hole, 
the photons may be accreted before they can diffuse out through the accre ting gas when th e 
accretion rate is high; they are trapped and advected into the black hole (iBegelman Ill979l ). 
Because the X-rays are most likely to be emitted very close to the black hole, more of them 
could be trapped than optical/UV photons, resulting in a steep aox- 



5.2. The Black Hole and Accretion Disk in PHL 1811 

It is generally thought that the source of the optical-UV bump in AGN broad band 
spectra is thermal emission from the accretion disk that powers the active nucleus. Despite 
this conviction, results from observations of the optical-UV properties in AGN do not conform 
to our expectations of accretion disks; specifically, the continuum does not have the expected 
slope, the dis k emission does not appe ar to be polarized, and there is no prominent Lyman 



edge feature (IKoratkar &: Blaes Ill999l ). PHL 1811 has a prominent big blue bump (Fig. 



in this section we explore the accretion disk explanation for this feature. 

The simplest model of an accretion disk is constructed by assuming that half of accretion 



energ y heats the disk, and the disk radiates locally like a black body (e.g.. lFrank. King &: Raine 



19921 ). This model has the advantage that it is easy to compute, but the disadvantage that is 
not thought to be physically realistic. As shown below, however, we can obtain some useful 
results by quantitatively comparing this model with our data. 



We compute the spectra emitted by the accretion disk between i? = 3.1 Rs and 1000 Rs 
in the frequency range log(z/) = 14-19 (Hz), for a range of black hole masses between 10^ 
and 10^^ M0, and a range of specific accretion rates, •m=Mj M^dd^ between 0.1 and 10, noting 
that this simple geometrically thin, optically thick model should break down at the higher 
specific accretion rates. We assume that the disk is observed face on (cosi = 0). To compare 
the accretion disk spectra with our data, we identify regions of the merged optical and UV 
spectrum that appear not to be dominated by emission lines, and obtain the average flux 
in these band^. We then compute the sum of the mean deviation of the log of vLy at the 
mean frequencies for each band from each model. The contours of the mean deviation are 



i°The bands we use are 1089-1102, 1317-1352, 1460-1482, 2229-2244, 3011-3038, 4014-4058, 4698-4768, 
and 5700-5750. 
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shown in the left panel of Fig. 9. We find that only a relatively small range of Mbh and 
rh match the data. The best fitting values are Mbh = 2.2 x IO^Mq and m= 0.9. The best 
fitting continuum spectrum is seen in the right panel of Fig. 9. 

What is the mass of the black hole in PHL 1811? We can estimate the mass using any 
number of the relationships betw een mass, H/3 velocity width, an d uLiy at 5100A currently 
available. We use the one given by lVestergaard fc Peterson I (120061 ). We measure ALa(5100A) 
to be 3.6 X 10^^ erg s~^ using a distance of 936.4 Mpc. We note that the flux from the merged 
spectrum appears to be consistent with the broad-band photometry. The mean flux density 
over a 980A band centered at observed 4400A is 1.0 x 10~^^ erg s~^ cm~^A~^. The photometry 
value, B = 14.4, corresponds to 1.15 x 10~^'*ergs~^ cm~^A~^, and so is consistent within 
15%. Given that the observations were not contemporaneous, this is quite good agreement. 
As will be discussed in P aper II, we measure the width of H/? to be 1943 km s~^. Using 
Vestergaard fc Peterson I (120061 ) Eq. 5, we find the black hole mass to be 1.8 x 10* M©, a 
factor of 12 smaller than the estimate based on the accretion disk model above. Thus, 
the disk is radiating at a rate significantly higher than the Eddington value to attain the 
luminosity in the optical-UV band pass that we observe. To illustrate this, we plot the 
spectrum for Mbh = 1-8 x 10* and ?fi=1.0 on Fig. 9. 

It is worth pausing to note that there are several reasons why this may not be an accurate 
estimate of the bla ck hole mass in very pecul i ar ob j ects such as PHL 1811. As is discussed in 
the introduction of I Vestergaard fc Peterson I (120061 ) and references therein, there are a num- 
ber of assumptions and suppositions that go into the construction of reverberation-based 
black hole masses. The primary assumption is that the line widths are virial. The evidence 
that this is true comes from a very few Sey fert 1.5 galaxies with excellent multiwavelength 
monitoring data (jPeterson fc Wandel Il2000l ). Narrow- line Seyferts Is like PHL 1811 are dif- 
ferent from these objects in that their H/? lines do not seem to vary, although variability has 
been observed in low- luminosity NLSls such as NGC 4051. The other evidence that the line 
widths are virial comes from the fact that in lower-luminosity reverberation-mapped AGN, 
the black hole mass estimated from the line widths and continuum luminosity is consistent 
with that estimated from the host galaxy stellar velocity dispersion and the Mbh^ct relation- 
ship. PHL 1811 is much different than these objects as it has a much higher luminosity. The 
uncertainty in the black hole mass in the reverberation-mappe d objects is given to be a factor 
of ~ 2.9 based on the scatter around the Mhh-cr relationship (lOnken et al. II2004I ). However, 
the deviation from the reverberation- mapping black hole mass in PHL 1811 may be larger 
for several reasons. First, it is possible that the broad line regions in Narrow-line Seyfert 



Is hav e a flattened configuration; this has been previously suggested by iMcLure fc Dunlop 
(I2OO2I ). A flatter distribution in NLSls would decrease the observed velocity width, and 



may also naturally yield the observed lower equivalent widths as a consequence of a smaller 
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covering fraction. This might imply a larger black hole mass in PHL 1811 than estimated 
by the regressions based on reverberation-mapped AGN. On the other hand, PHL 1811 has 
a soft spectral energy distribution. This means that it produces a weaker ionizing contin- 
uum compared wit h other quasars w i th sim ilar optical luminosity. Inverting the argument 
originally made by IWandel fc Boiler I (119981 ) (applied to NLSls that were inferred to have a 
stronger ionizing continuum because of their prominent soft excess), this implies that the H/3 
lines should be produced closer to the continuum source than in other NLSls, which would 
require the black hole mass to be smaller. The point is that extrapolating the relationships 
obtained from the reverberation-mapped AGN may be particularly dangerous for extreme 
objects such as PHL 1811. 

Regardless, assuming the black hole mass estimated above, we infer that the disk must 
be radiating at a rate very much higher than the Eddington value to attain the observed 
optical-UV luminosity. We n ote we are not the first to in fer super-Eddington radiation in 
NLSls; this was also found by lCoUin fc Kawaguchi I (120041 ) using a different approach. They 
estimate the black hole masses from the H/? line widths, and then compare the observed 
bolometric luminosity to the predicted, whereas we instead compare the luminosity and 
shape of the accretion disk to the observed continuum. 

What is the bolometric luminosity of PHL 1811? We integrate over the inferred broad- 
band continuum to determine this. For wavelengths shorter than one micron, we use the 
spectral energy distribution inferred from the observed optical, UV and X-ray data. It is 
comprised of piecewise power laws, suitable for use in Cloudy, a nd ca n be seen in Fig. 12 
of Paper II. Longward of one micron, we use the lElvis et al. I (119941 ) average continuum 
spectrum, since it is seen to match the 2MASS and IRAS photometry rather well (Fig. 
8). The result is 3.7 x lO^^ergs^^ It is 10.2 times ALa(5100), and so U^i/ \Lx{^lQQ) is 
very close to 9, the bolometric correction factor ofte n used for AGN (and assumed in the 
mass- luminosity relationship in iPeterson et al. 1120041 ). 



The Eddington luminosity for a 1.8 x 10® M© black hole is 2.25 x 10^^ erg s~^. This implies 
that PHL 1811, overall, is radiating at about 1.6 times the Eddington luminosity. But as the 
accretion disk fits show, the optical-UV region is radiating at a very super-Eddington rate. 

The sum-of-black bodies spectrum is not a physically realistic model for an AGN ac- 
cretion disk, although since it is the sum of black bodies, it will be the overall brightest disk 
for a homogeneous, isotropically emitting disk. It does not take into account vertical disk 
structure, Comptonization, or other effects that are expected to influence the disk spectrum. 
What do we expect to see from a realistic accretion disk spectrum in comparison with the 
sum of black bodies model? 
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Shimura fc Takahara I (jl993l . Il995l ) considered the effect of vertical structure and radia- 
tive transfer in tlie accretion disk, and applied their results to AGN and black hole candidates. 
They find that the emerging spectrum is significantly harder than the sum of black bodies 
disk, due to electron scattering and the vertical temperature gradient. This means that for a 
given accretion rate and black hole mass, their computed disk spectrum has proportionally 
more soft X-rays than the sum-of-black-bodies disk, but at the expense of the UV. In order 
for the flux of such a disk spectrum to match the UV flux from PHL 1811, an even larger 
black hole would be necessary. In fact, this shou ld be generically true for any disk mode l 



that radiates at less than the Eddington rate (e.g.. lKawaguchi. Shimura fc Mineshige II2OOII ). 



More recently, mode ls for disks with super-Ed dington ern i ssion have been proposed 



by several authors (e.g., iBegelman 1 12002| ). In the iBegelman I (120021 ) model, the disk is 



inhomogeneous, and radiation leaks out through optically thin channels. We infer that 
PHL 1811 is radiating at a super-Eddington rate in the optical and UV. Perhaps this means 
that the spectrum of a super-Eddington disk should be similar to that of a typical quasar 
longward of ~ 1550A, based on Fig. 8 of Paper II. One might expect the super-Eddington 
emission to be stronger closer to the black hole, where the shorter wavelength emission 
originates. Perhaps this is the origin of the unusual rising continuum shortward of ~ 1400A 
shown in Fig. 8 of Paper II. 



5.3. The Host Galaxy of PHL 1811 



Magorrian et al. 



(1l998l ) discovered that black hole masses are correlated with the lu- 
minosity of the host bulge or elliptical galaxy. Later it was found that a tighter correlation 
between th e velocity dispersion of the bulge of a galaxy and the mass of the nuclear black 
hole exists (IFerrarese &: Merritt II2OOOI : iGebhardt et al. Il2000l ) . The host galaxies of luminous 
active galaxies, i.e., quasars, are almost exclusively ellipticals, and their basic properties are 



the same as elliptical galaxies without quasars (e.g.. lDunlop et al. II2003I ). Furthermore, they 



obey the relationship b etween black hole mass and spheroid mass observed in nearby galax- 



ies (e.g.. iDunlop Il2004h . PHL 1811 has a rather large black hole (~1.8x10®Mq; §5.2) and 
therefore should have a large elliptical galaxy as a host. An object that should in principle be 
quite similar is t he quasar 3C 273. Its re dshift is 0.158 and it s black hole mass is inferred to 
be 8.8 X 10^ M0 JPeterson et al. Wim^ . iMartel et al~l J^OoJ present the ACS coronagraph 
data from this object. They find a large elliptical galaxy with an inner region about 17 kpc 
in diameter (6.5"), and extended diffuse emission out to 6-12" from the active nucleus. A 
de Vaucouleurs fit to the V & I profiles yield effective radii of ~ 2.2" and ~ 2.6". 



Jenkins et al 



((20051) show an HST ACS WFP image of PHL 1811. Observations of the 
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host galaxies of nearby quasars are best done using the coronograph; however, there was not 
sufficient time in that program to make the coronagraph image with an accompanying PSF 
cahbration exposure. They instead made a 520 second observation directly using the F625W 
(Sloan r*) filter, and corrected for the QSO PSF using a nearby star with colors similar to 
those of PHL 1811. The resulting image was surprising; it does not show a large elliptical 
host like 3C 273; rather there are structures that look like spiral arms on either size of the 
QSO. These features are not an artifact of the PSF subtraction as they can even be see n 
in the image before subtraction, and for other reasons discussed by iJenkins et al. I (120051 ). 
Thus, PHL 1811 appears to have a spiral host, and appears to lack a large bulge that might 
be expected given its luminosity. 

The 3C 273 observations were 9-10 times longer than the PHL 1811; perhaps the 
PHL 1811 image is not sufficiently well exposed to observe the bulge? We can estimate 
the expected properties of a typical host galaxy, by using our black hole mass estimate to 
estimate the mass of the spheroid, then determining the expected observational properties 
of an elliptical galaxy having a spheroid of that mass. 

In §5.2 we showed that the luminosity of the optical-UV spectrum is consistent with 
a black hole of 2.2 x lO^ MfD. However, the black hole m ass based on AGN reverberation 
mapping results given by IVestergaard &: Peterson I (120061 ). we obtain a much smaller black 
hole mass of 1.8 x 10^ Mq. Using the r elationship bet ween the black hole mass and the mass 
of the spheroid {Mbh = 0.0012M,p,,; IPunlop I Eool yields Msph = 1-5 x W^Mq for the 
estimated Mbh = 1-8 x 10^ M©. 

It has been found that the host galaxies of quasars are indistinguishable from elliptical 
galaxies. Therefore, we can use the properties of elliptical galaxies to estimate the size and 
brightness of the putative host. The mass of the sphero id is related to the effective r adius by 
re = 2.6(M,p/,/10")3/5kpc for Hq = 70kms-^Mpc-i JLoewenstein fc White lEoogl ). Thus, 
for PHL 1811, we expect an effective radius of 3.3 kpc. At the distance of PHL 1811, one 
arc second corresponds to 3.2 kpc, so rg = 3.3 kpc corresponds to 1". The image presented 
by iJenkins et al. I (120051 ) is not reliable within about 1.5" due to the correction for the QSO 
PSF, but should be fine at ~ 2re. 

It is well known that the surface brightness of elliptical galaxies is correlated with the 
effective radius (e.g., 



Diorgovski fc Davis 1119871 ). There is, however, considerable scatter. 



although iMcLure fc Dunlop I (120021 ) argue that the rms scatter for AGN hosts is only 0.18 
dex. We us e the results from a r ecent analysis of elliptical galaxies in the Sloan Digital 
Sky Survey (IBernardi et al. II2003I ). For the effective radius of 3.3 kpc, their Fig. 10 shows 
that the mean surface brightness is ~ 20 magarcsec"^, and almost all objects lie within 19- 
21 magarcsec"^ in the r* filter. We need to account for cosmological dimming, which will 
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lower the observed surface brightness by a factor (1 + z)^ = 2.02. Then, the mean surface 
brightness is 20.75 mag arcsec"^, with a conservative range of 21.75-19.75 mag arcsec"^. 

Since the image is not rehable at 1 re due to psf distortion, we need to estimate the 
surface brightness at 2 rg. Assuming a de Vau couleurs profile, we can use this average to find 
the expected surface brightness at 2 (e.g.. iBinney fc Merrifield Ill998l ). We find then the 
mean surface brightness to be 22.25 mag arcsec"^, with a range of 23.25-21.25 mag arcsec"^. 

Next, we use the exposure time calculator for the HST ACS to determine the ex- 
pected signal-to-noise ratio of diffuse emission obtainable in a 520 second exposure using 
the F625W filter. We use the built-in spectrum of an elliptical galaxy, and include the effect 
of Galactic extinction by E{B — V)= 0.046. We find that the mean surface brightness of 
22.25 mag arcsec"^ yields a signal-to-noise ratio in a 2x2 pixel box of 5.4, with a range of 
2.3-11.7 for the surface brightness range. Thus, if PHL 1811 had a norm al elliptical galaxy 
host, it may have been detected at 2 re = 2" in the lJenkins et al. I (120051 ) image, assuming it 
is on the brighter side of the range of elliptical galaxies. 

Finally, it has been found that NLSls frequently have a much larger host galaxy bulge 
than expected based on their black hole mass estira ated from th eir emission- l ine w idths 
fe.g.- lGrupe fc Mathur Ibooi iMathur fc Grupe Il2005h . In addition. [Rvan et al. I J2006h find 



from near-IR imaging of NLSls that the bulge is about an order of magnitude larger than 
would be inferred using the Mbh ^ cr relationshiop and the black hole masses estimated from 
reverberation mapping. If PHL 1811 had followed this trend, the host galaxy bulge should 
have been very easy to see. 



6. Summary 

In this paper, we report the results and analysis of two Chandra observations that were 
coordinated with an HST observation, an XMM-Newton observation, two Swift ToOs, and 
MDM optical photometry of the unusually luminous, nearby narrow-line quasar PHL 1811. 
Here we summarize the primary results of the paper. 

• The two lOks Chandra observations, made twelve days apart, reveal a weak X-ray 
source at the position of PHL 1811. A factor of four variability was observed between 
the two observations. The X-ray spectrum is steep, with photon index F ~ 2-2.6, 
typical of the power law indices observed in ASCA spectra from NLSls. There is 
marginal evidence of a soft excess in the brighter spectrum, and marginal evidence for 
spectral variability in which the soft excess became more prominent when the object 
was brighter. No evidence for absorption was found. 
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• The ~ 30 ks XMM-Newton observation revealed no evidence for X-ray variability, and 
a flux consistent with the fainter of the two Chandra spectra. The X-ray spectrum 
is steep with F = 2.3 ± 0.1, again typical of the power law indices observed from 
NLSls by ASCA. There is no evidence for absorption, with an upper limit on intrinsic 
absorption of 8.7 x 10^° cm~^. The Optical Monitor data revealed no evidence for short 
time-scale UV variability. Two Swift observations found the flux to be consistent with 
the brighter of the two Chandra spectra. Overall, a factor of ~ 5 variation in X-ray 
flux was observed among the flve X-ray observations presented in this paper. The UV 
photometry, obtained by XMM-Newton and Swift, is consistent with the i/5T spectrum 
discussed in Paper II; thus no evidence for UV variability has been found. In addition, 
MDM photometry observations revealed no evidence for optical variability over three 
nights. 

• We compare the XMM-Newton PN spectrum of PHL 1811 with those from two repre- 
sentative NLSls: Ton S180, which has been previously modelled using a dual Comp- 
tonization model, and IH 0707—495, which has been previously modeled using partial 
covering or a reflection model. Since the comparison NLSls are much brighter than 
PHL 1811, we extracted spectra from sufficiently short segments that the number of 
photons in the PN spectra are approximately the same as in the PHL 1811 spectrum. 
Fitted with a power-law model, the comparison spectra are significantly steeper than 
that of PHL 1811, and complexity is detected in the spectrum from IH 0707—495. 
PHL 181 1's X-ray spectrum seems to be different from those of other NLSls. Since 
the slope is consistent with the power- law observed in ASCA spectra of NLSls, the 
simplest explanation for the X-ray spectrum in PHL 1811 is that it is powered by 
Compton upscattering of soft photons by energetic electrons, the typical power-law 
X-ray emission mechanism in quasars. 



The i ?epiPo6'y4X observation found PHL 1811 to be significant X-ray weak (ILeighly et al 



20011 ). At that time, we presented three possible explanations for its X-ray weakness: 
1.) PHL 1811 is a BALQSO; 2.) PHL 1811 is temporarily in a low fiux state; 3.) 
PHL 1811 is intrinsically X-ray weak. The spectral energy distribution, constructed 
using the C/iandra spectra and the /fST spectrum taken two days before the first Chan- 
dra observation, confirms that PHL 1811 is X-ray weak, with aox between —2.2 and 
—2.4. The XMM-Newton and Swift observations are consistent with these values. The 
aox for a typical quasar of similar UV luminosity is —1.6, and compared with a sample 
of optically-selected quasars of similar UV luminosity, the X-ray fiux of PHL 1811 is a 
factor 13-450 times weaker. The lack of absorption in the X-ray spectra, and the lack 
of absorption lines in the UV spectrum (Paper II) conclusively rules out the possibility 
that PHL 1811 is a BALQSO, with a weak, absorbed X-ray spectrum characteristic of 
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that class of object. PHL 1811 has now been observed seven times in X-rays between 
1990 (during the ROSAT All Sky Survey) and the most recent Swift observation 2006 
May 14. While it varies by about a factor of 5 among the latter 5 observations, ruling 
out a scattering origin for the X-ray emission from an extended electron-scattering 
mirror, it is always significantly weaker than other quasars of similar optical lumi- 
nosity. While we can never rule out the possibility that PHL 1811 is coincidentally 
observed always when it is in a transient low state, the plausibility of this hypothesis 
is decreasing. We conclude that PHL 1811 is intrinsically X-ray weak. 

We discuss the possible origins of the X-ray weakness. It may be a consequence of the 
high accretion rate inferred from the narrow lines and high luminosity. The corona may 
be smaller, and with a larger fraction of the accretion power going into the optically- 
thick, geometrically-thin accretion disk emitting the optical and UV continuum. Alter- 
natively, the corona may be quenched by the strong optical and UV emission, and may 
be unable to Compton-upscatter those photons to X-ray energies. Another possibility 
is that high accretion rate traps the X-ray photons and advects them into the black 
hole. 

We compare the spectral energy distribution with a simple sum-of-black bodies accre- 
tion disk model. We find that the UV spectrum is consistent with a black hole mass 
of 2.2 X 10^ M0 and a specific accretion rate of m= 0.9. Using the H/3 FWHM, the 
luminosity at 5100A and an equation derived from results of reverberation mapped 



AGN given by IVestergaard fc Peterson I (120061 ). we obtain a black hole mass estimate 
of 1.8 X 10^ Mq, a factor of 12 below the value obtained from the accretion disk spec- 
trum. We caution that we do not have a direct measurement of the black hole mass for 
this object. We point out that while the sum-of-blackbodies disk is not physically real- 
istic, it is the brightest homogeneous, isotropically-emitting disk, and most other, more 
realistic, disks will require an even higher radiation rate relative t o Eddington. The 



exception may be inhomogeneous disks recently proposed by e.g., iBegelman I (120021 ) 



in which radiation leaks out through optically-thin channels. Based on t he observed 



broad-band continuum spectrum, and extrapolating to the IR using the lElvis et al. 



( I1994J ) SED, we estimate the bolometric luminosity to be 3.7 x lO'^^ergs ^. Thus, over- 



all, PHL 1811 is radiating at 1.6 times the Eddington rate, with the optic al and UV well 



above Eddington. We also discuss the HST ACS image presented by IJenkins et al. 



( 120051 ) that reveals a spiral galaxy at the position of PHL 1811. Most quasars have 
elliptical hosts; we demonstrate that if PHL 1811 had the large elliptical host typical 
of average quasars with 1.8 x IO^Mbh black holes, it might have been seen (i.e., the 
estimated SNR is ~ 5) in that image. A spiral host galaxy appears to be yet another 
unusual feature of PHL 1811. 
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Table 1. Observing log 



Observatory andr 


Date 


Exposure 


Bandpass or 


Extraction Radius 




Instrument 


(seconds) 


Effective Wavelength 


(aJTcseconds) 


Chandra ACIS S3 


2001 Dec 5 


9377 


0.3-10 keV 


3.94 


Chandra ACIS S3 


2001 Dec 17 


9839 


0.3-10 keV 


3.94 


XMM-Newton PN 


2004 Nov 01 


27472 


0.3-12 keV 


27 


XMM-Newton MOSl 




32126 


0.5-10 kcV 


23 


XMM-Newton MOS2 




32164 


0.5-10 kcV 


23 


XMM-Newton OM (UVM2) 




25800 


2310A 


12 


XRT (PC mode) 


2005 Oct 22 


2462 


0.3-10 keV 


23.4 


5OT/t UVOT (UVWl) 




787 


2600A 


12 


ftLift UVOT {UVM2) 




844 


2200A 


12 


Swift UVOT (UVW2) 




844 


1930A 


12 


Swijt XRT (PC mode) 


2006 May 12 


1600 


0.3-10 keV 


23.4 


Sm/i UVOT (V) 




130 


5460A 


6 


Swift UVOT (B) 




130 


4350A 


6 


Swift UVOT (U) 




130 


3450A 


6 


5m/C UVOT (UVWl) 




259 


2600A 


12 


Swift UVOT (UVM2) 




372 


2200A 


12 


UVOT (UVW2) 




526 


1930A 


12 


MDM McGraw-Hill 1.3 + I + Tcmplcton 


2004 Oct 14-16 


bAl^ 


8050A 


5.8-10.0 


MDM McGraw-Hill 1.3 -|- V -|- Templeton 




5,9,5^ 


5380A 


5.7-10.1 


MDM McGraw-Hill 1.3 -h B -|- Templeton 




5,5,2^ 


4350A 


6.2-10.0 


MDM McGraw-Hill 1.3 -t- U -(- Templeton 




5,5,2^ 


3640A 


6.8-10.2 



''Number of frames on each of the three days. 
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Table 2. Chandra Spectral Fitting Results 



Parameter 2001-12-05 2001-12-17 



1. Power law model, spectra fit separately. 



Photon Index 
Normalization*^ 

Flux (0.3-5 keV; ergcm~^s~^) 
Luminosity (0.3-5 keV; ergs~^) 



-0.37 
-0.36 



2.01: _ 

1.17 ±0.02 X 10-5 
5.2 X 10-^^ 
5.6 X 10^2 
1.96/2 



^•'-'0-0.18 

5.10 ±0.45 X 10-5 

2.3 X 10-^^ 

2.4 X 10^3 
13.2/15 



2. Power law ± blackbody'^ 
Photon Index 
PL Normalization^ 
Blackbody temperature (keV) 
BB Normalization'^ 
xVd.o.f 



2.22 ±0.34 

AAtlf X 10-5 

096+°-°^° 
u.uyu_Q 

8.36/13 



3. Power law ± blackbody*^ , joint fit 

Photon Index 2.12 ±0.25^^ 

PL Normalization^ 1.17 ± 0.22 x 10-^ 
Blackbody temperature (keV) 0.10l[^:[jl^<^ 
BB Normalization^ 0^ 
xVd.o.f 10.8/16^^ 



2.12 ±0.25^^ 
4.14+°;^^ X 10-5 

-l-0.035d 
-0.042 

i.5ig:L X 10-' 



0.10 



10.8/16'^ 



^In units of photons keV"^ cm^^s"^ at 1 keV in the observed frame. 

'^The blackbody temperature is given in the rest frame. 

^In units of L-^g/Dl^, where L39 is the source luminosity in units of 
10^^ erg s-^, and Diq is the distance to the object in units of 10 kpc. 

'^Parameter fit jointly. 

^ Fixed parameter. 
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Table 3. XMM-Newton Spectral Fitting Results 



Parameter Measurement 



Photon Index 


9 90+0.12 


PN Normalization'^'* 


1.16 ±0.98 


MO SI Normalization Offset'' 


1 9P;+0.21 


M0S2 Normalization Offset'' 


1 09+0.22 


Flux (0.3-5 keV; ergcm~^s~^) 


5.1 X 10-1^ 


Luminosity (0.3-5 keV; ergs^^) 


5.4 X 10^2 


xVd.o.f. 


74.7/79 



^The spectral model was Ce^H{Gai)a{E)^-r ^ 
where all parameters were tied together except for 
the constant C, which was fixed to 1 for the PN 
and allowed to vary for the MOSl and M0S2. Un- 
expected optical loading plausibly causes the MOS 
spectra to have a higher normalization; therefore 
we quote flux and luminosity values from the PN 
data only. 

^In units of photons keV~^ cm~^ s~^ at 1 keV in 
the observed frame. 
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Table 4. Swift UVOT Results 







Flux Density (10" 


i4gj,gg 1cm ^ A 




Observation 


UVW2 


UVM2 UVWl 


U B 


V 


2005 October 22 


3.75 ±0.03 


2.75 ±0.02 2.45 ±0.02 






2006 May 12 


3.66 ±0.03 


2.71 ±0.03 2.47 ±0.03 


1.89 ±0.04 0.99 ±0.02 


0.61 ±0.02 



'^Observed fluxes; uncorrected for Galactic reddening. 
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Fig. 1. — The Chandra image {left) and the MDM optical i2-band image (right) show that 
the X-ray source is securely identified as PHL 1811. 
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Fig. 2. — Top: Unfolded best fitting power-law models for the spectra from the two Chandra 
observations. Bottom: Residuals from a power-law model fit to the spectrum from the second 
observation have a convex shape suggesting the presence of a soft excess. Inset: The ratio of 
second spectrum to the first spectrum shows that the spectral variability is predominately 
confined to the softest energies, suggesting the emergence of a soft excess component when 
the object is bright. 
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Fig. 3. — XMM-Newton background-subtracted lightcurvc composed of the 0.5-5 keV pho- 
tons from the MOSl, M0S2 and PN detectors. The lower hght curve shows the estimated 
background in the source extraction region. There is no strong evidence for variabihty over 
the 31 ks observation. 
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Fig. 4. — The XMM-Newton MOSl, M0S2 and PN spectra (open squares, open triangles, 
and filled circles, respectively) are well fit using a power law and Galactic absorption. The 
upper limit on intrinsic X-ray absorption is 8.7 x 10^°cm~^. 
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Fig. 5. — A comparison of the XMM-Newton PN spectrum from PHL 1811 with spectra 
from other NLSls from segments of data short enough for the statistics to be similar. In 
each panel, the ratio of the data to a power law plus Galactic absorption model is plotted. 
The resulting power law index and xl given in each panel. Top: PHL 1811 is statistically 
well described by a power law model. The p hoton index is ty pical of the power law indices 



observed from ASCA observations of NLSl s (ILeighlv 



1999bf ). Middle: Ton S180, an NLSl 



classified as having a "simple" spectrum by I G alio I (120061 ). is statistically well described by a 
power law; how 'ever, the index is stee p, suggesting that a longer exposure would reveal a hard 
tail (as it does, Vaughan et al. 2002h . Bottom: IH 0707—495, an NLSl classified as having a 



Gallo 



(120061 ) ■ shows a ver y steep spectrum an d significant residuals 



"complex" spectrum by 
confirming the necessity of a complex model (e.g.. lGallo et al. Il2004bl : iTanaka et al. 1120041 ). 
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Fig. 6. — Long-term X-ray flux variability of PHL 1811 in terms of uFi, at 2 keV in the 
rest frame. The uncertainties for the spectroscopic data {Chandra and XMM-Newton) are 
propagated 1-a errors in the power law normalization and index. For the detection data 
[Swift) the uncertainties are proportional to the count rate error. Note the logarithmic flux 
axis. Also shown is the predicted flux for aox = —1.6 based on the UV flux from the HST 
observations. 
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Fig. 7. — The XMM-Newton OM and Swift UVOT photometry points are overlaid on the 
observed merged optical and HST UV spectrum observed 2001 Dec 03 and discussed in 
Paper II. No strong evidence for UV variability is seen among the observations which span 
4.5 years. 
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Fig. 8. — The spectral energy distribution of PHL 1811, plotted as a function of the rest- 
frame frequency. Contours from each of the two Chandra observations are shown, generated 
by successively setting each parameter to its Ax^ = 2.71 value and computing the model, 
then determining the maximum and minimum of all of the models. The dashed line shows the 
ex pected 2 keV flux for a n average quasar of this luminosity, based on the regres sion presented 
bv IWilkes et al. I (119941 ) , while the do tted line shows th e range observed by IWilkes et al. 
( I1994J ). The average quasar SED from lElvis et al. I (119941 ) . scaled to the 1 micron inflection 
is also shown. 



-45 - 




og Black Hole Mass 




log Frequency (Hz) 



Fig. 9. — Left: Contours of equal deviation of UV flux points from sum-of-black bodies 
accretion-disk models, as a function of input black hole mass (in units of solar masses), and 
the speciflc accretion rate, m. See text for further information. The contour interval is 0.05 
in units of loguLy. The minimum is located at Mbh = 2.2 x 10^ Mq and m=0.9. Right: 
The UV continuum flux points and the X-ray spectra are shown as solid dots and solid 
lines respectively. The best fltting model is shown by the dashed line. For comparison, the 
spectrum for a Mbh — 1-8 x 10^ Mq black hole with m^l.O is shown by the dotted hne. 



